Ambipolar transition voltage spectroscopy: analytical results and experimental 
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This work emphasizes that the transition voltages Vt± for both bias polarities (V<0) should be 
used to properly determine the energy offset eo of the molecular orbital closest to electrodes' Fermi 
level and the bias asymmetry 7 in molecular junctions. Accurate analytical formulas are deduced 
to estimate £0 and 7 solely in terms of Vt±. These estimates are validated against experiments, by 
showing that full experimental J-V-curves measured by Beebe et al [Phys. Rev. Lett. 97, 026801 
(2006)] and Tan et al [Appl. Phsy. Lett. 96, 013110 (2010)] for both bias polarities can be excellently 
reproduced. 

PACS numbers: 85.65.+h, 85.35. Gv, 73.63.-b 
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I. INTRODUCTION 

Transition voltage (Vt) spectroscopy (TVS) has been 
recently proposed^ to get insight into the energy offset e 
between the metal Fermi energy and the closest molec- 
ular orbital (HOMO or LUMO), which has long been 
recognized to be a key quantity for the charge trans- 
port in molecular devices 2 -. Due to its simplicity, it soon 
became very popular for interpreting molecular trans- 
port measurements.— ~— Within the initial proposal re- 
lying upon a tunneling barrier picture^ Vt, the mini- 
mum of the Fowler-Nordheim (FN) plot \n(I /V 2 ) versus 
1/V determined from I-V measurements, has been asso- 
ciated to the point where the barrier tilted by the applied 
bias change from trapezoidal to triangular. This yields 
eVf = £0^ Later, TVS was interpreted within a coherent 
transport model based on a single levelJ^ In the (real- 
istic) cases where the energy offset is sufficiently larger 
than the level broadening due to the couplings to elec- 
trodes, the relationship deduced within the latter model 
(eVt = 1.15eo for a symmetrical orbital location between 
electrodes^) turned out to be not much different from 
the original "barrier shape" conjecture. With certain lim- 
itations, ab initio studie o 17 i 18 give microscopic support to 
the single-level model. A significant aspect in the TVS 
analysis^ is to properly account for the potential pro- 
file asymmetry. In the presence of this asymmetry, the 
FN-plots also become asymmetric, a fact which reflects 
itself in different magnitudes of the transition voltages 
Vt+ 7^ — Vt- for both bias polarities (V<0), as pointed 
out recently^ It is the main aim of this paper to consider 
the ambipolar TVS in detail. 

The remaining part of this paper is organized as fol- 
lows. The general theoretical framework will be pre- 
sented in Sect. [TTJ In Sect. IIII1 accurate analytical for- 
mulas will be given enabling one to directly extract the 
quantities of physical interest from the transition voltages 
measured for positive and negative biases. The accu- 
racy and the usefulness of these theoretical formulas will 
be illustrated in Sect. IIV1 by showing that applications 
to experimental /-V-data yield an excellent agreement. 
Conclusions will be presented in Sect. [V] 



II. THEORETICAL FRAMEWORK 

An applied bias V affects the energy offset eo = 
£o(^OIv=o — > £o(V). Most easily, as currently done in 
electrochemistry^ - — this effect can be accounted for by 
means of a voltage division factor— 7, which specifies the 
bias at the location of the orbital's "center of gravity" 



eo(V) = s +jeV. 



(1) 



With a potential origin as in Fig. [T] 7 ranges from —1/2 
to +1/2. Eq. ([1} corresponds to a potential that is flat 
across the molecule and entirely drops at contacts. The 
intcrfacial potential drops are 

SV S = V/2-[e (V)~e ]/e=(l/2- 1 )V 

SV t = [e (V)-e ]/e + V/2 = (l/2 + 7)V. (2) 

A positive (negative) 7-value corresponds to a larger 
(smaller) potential drop at the positive electrode, or, al- 
ternatively, to a molecular orbital energy shifted upward 
(downward) by a positive bias. 

In the wide-band limit, wherein the transmission is 
Lorentzian, the current through a single level (Newns- 
Anderson mo del) ca n be expressed analytically as (see, 

e. g ., Refs. mmm 



2e T 2 
I = N tltA 

h r„ 



A 



A_ 



arctan — arctan — — I . (3) 



Here, N is the (effective) number of molecules contribut- 
ing to the c urrent , A± = e (V)±eV/2. T a = (T s + T t ) /2 
and T g = ^/r s r t , T s _ t being the level broadenings due to 
molecule-electrode couplings. In usual cases of interest, 
r a |e | and voltages not much higher than V t r^ the ar- 
guments of the inverse trigonometric functions of Eq. ^ 
are large, and one can approximate 
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The comparison with the results based on the exact 
Eq. ([3]) shows that Eq. (j4]) very accurately describes the 
FN-transition; see Fig. [21 
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FIG. 1: Schematic representation of a HOMO-mediated con- 
duction (the HOMO much closer to the Fermi level than the 
LUMO, Eh "C si) in a setup characterized by asymmetric elec- 
trodes [asymmetric voltage division y(> 0) and asymmetric 
molecule-electrode couplings F s and Ft]. See the main text 
for details. 
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FIG. 2: Fowler-Nordheim (FN) plots for a single-level coher- 
ent transport obtained by means of the exact Eq. ([3]) and 
approximate Eq. (T3J for symmetric (7 = 0) and asymmetric 
(7 7^ 0) bias profiles. They show that for level broadenings F a 
sufficiently smaller than its energy offset (eo = 1 eV), the FN- 
transition is accurately described by Eq. Q. See the main 
text for details. 



III. USEFUL ANALYTICAL RESULTS FOR 
TRANSITION VOLTAGE SPECTROSCOPY 



Eq. ((3]) can be used to deduce simple analytical ex- 
pressions of the relevant quantities within the (realistic) 
assumption of a level broadening sufficiently smaller than 
the energy offset, which are exact to O (V a /eo) 2 . By im- 
posing d\n(I/V 2 )/d(l/V) = 0, one obtains the transi- 



Xti = eo/Vti = -2 7 + ^7 2 + 3/4, 
Xti = e /V t2 = -2 7 - V7 2 + 3/4. 



(5) 



Because —1/2 < 7 < 1/2, it is easy to show that 
Xti > and Xt2 < 0. Therefore, the signs of V t \ and 
Vt2 are opposite. Denoting by V t +(> 0) and Vt-(< 0) 
the transition voltage for positive and negative polar- 
ities, Vt+ = Vti and V t - = V t 2 for LUMO-mediated 
transport (eo > 0), while for HOMO-mediated transport 
(s = Sh < 0) V t + = V t2 and V t - = V n . In the HOMO 
case, V t + < |V r t _| for 7 > 0, whereas V t+ > \V t -\ for 
7 < 0. For 7 = 0, the result \V t ±/e \ = 2/V3 = 1.15 for 
symmetric case^ is recovered. 

Concerning the signs in general, it is worth noting that, 
according to Eq. (J]), a redefinition of the bias polarity 
(V — > —V) yields a sign change in the division poten- 
tial factor (7 — ^ —7). Therefore, the discussion can be 
restricted to the range, e. g., < 7 < 1/2. 

To illustrate the accuracy of the transition voltages V t ± 
expressed by the above analytical formulas, a compari- 
son with the transition voltages deduced from the exact 
Eq. ([3]) is presented in Fig. [3] As visible there, the ex- 
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FIG. 3: The reduced transition voltages Vt±/eh plotted versus 
the inverse of the relative HOMO-level (eo = — £h) broadening 
Fa/£h for given potential division factors (7 = 0,0.1) deduced 
exactly [Eq. ©, thick lines] and approximately [Eq. @, thin 
horizontal lines]. Notice the rapid saturation for sufficiently 
small level broadenings r a . See the main text for details. 



act results are rapidly approached for small T a /sa — > 0. 
Graphical results obtained by means of Eq. are pre- 
sented in Fig. |4] for values < 7 < 1/2. In the opposite 
case ( — 1/2<7<0) they can be deduced by symmetry 
V t± (- 7 ) = ±V tT ( 1 ). 

An important consequence of Eq. is that both the 
voltage division factor 7 and the energy offset eg can be 
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FIG. 4: (a) Results for the transition voltages Vt± as a func- 
tion of potential division factor 7 deduced from Eq. J5). (b) 
Ratio between the energy offset £0 = — £h < and the tran- 
sition voltages Vt± of both bias polarities as a function of the 
transition voltage asymmetry Vt-/Vt+. Each point of this 
curve corresponds to a given bias asymmetry 7, and a few 
values important for the discussion in the main text are in- 
dicated: 7 = 0.25 -»■ \Vt, max /Vi, min \ = 3.5, 7 = 0.067 -> 
|e /VU<«l = 0.40, |V t _| = |e |, and \Vt, mi n/Vt, min \ = 1.36. 
See the main text for details. 

TVS's proof of value for molecular electronics is the 
fact that the FN-minimum occurs at voltages below the 
values corresponding to resonant tunneling where the dif- 
ferential conductance is maximum. This is important be- 
cause, with seldom exceptions ; 14 ' 27 molecular junctions 
cannot withstand such high voltages. Still, as already 



noted) 4 -^ it is only a small range V > Vt (if at al&) 
that can be sampled in experiments. The situation can 
be further improved by using the minimum V = Vt±(n) 
of a generalized FN-plot ln(I/V K ) vs. l/V (1< K < 2)M 
General analytical expressions valid for arbitrary k can 
also be deduced 
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The assignment Vti^K) — > Vt±(«) is the same as dis- 
cussed above for n = 2. The usage of k- values smaller 
than 2 results in lower transition voltages V t ± (n),— which 
can easier be sampled experimentally. On the other side, 
being smaller they are more affected by the relatively 
large (> 0.1 V 12 ' 14 ) experimental errors. 

From a pragmatic standpoint, Eqs. ([5]) and (J7J as well 
as the more general Eqs. © and (fTU)) represent the core 
of this paper: they enable one to determine the quantities 
of physical interest £0 and 7 from the transition voltages 
Vt± measured for both bias polarities. 

IV. DISCUSSION OF EXPERIMENTAL DATA 

The above analytical results hold in general for the 
transport mediated by a single level whose energy offset is 
sufficiently larger than its electrode-induced broadening. 

I will next apply these results to the experimental data 
of Refs. Q] and Q3, obtained for the HOMO-mediated 
transport through molecular junctions in CP-AFM (con- 
ducting probe-atomic force microscopy and CW 
(crossed- wire) 1 - setups. For a concrete comparison with 
experiment, it is obviously important to correctly assign, 
out of the two transition voltages measured for opposite 
bias polarities in experiment, which is Vt+ and which is 
Vt-. Therefore, before entering into details, I note that 
the positive and negative biases have been chosen as the 
ones utilized in the experiments of Refs. [1IIT2I discussed 
below. 

Table |T] collects the experimental data for Vt± from 
Refs. [TJ and [l2|. They have been used to compute the 
values of en and 7 from Eqs. © and ([7]), which are also 
given in Table HI 

Further, I will use these values of £/, and 7 to re- 
produce the available I-V-A&ta, for anthracenethiol- and 
terphenylthiol-based junctions measured in Refs. [l| and 
[13^), respectively. To this aim, I will employ Eq. (QJ, 
which represents a very good approximation of the exact 
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Eq. ([3]) for biases not too much larger than the tran- 
sition voltages. The prefactor NT 2 g in Eq. ^ can be 
determined from the experimental linear conductance 
(N^Tg = 0.051 eV and iV 1 / 2 ^ = 0.124 eV, respec- 
tively). The theoretical curves obtained in this manner 
are plotted against the experimental ones in Figs.EK and 
b. The agreement is excellent, and this demonstrates the 
remarkable accuracy of the present approach based on 
Eqs. flU and 

Table U shows that for the molecular junctions of 
Refs. [H and [l2l the energy offset is, in fact, very close to 
the estimate of the barrier-shape conjecture if the transi- 
tion voltage for positive biases is employed, eVt+ = Eh- 6 
As visible in Fig. HJd, this does not hold in general but 
only for potential division factors close to 7 = 0.067. 
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FABLE I: Experimental values of Vt± for anth(racene)- and 
terphenyl (TP)-based molecules and platforms from Refs. [l] 
(a) and[l2 (b), and values of Eh and 7 calculated via Eqs. ([6]) 
and |7|. See the main text for details. 

Noteworthy, all the cases presented in Table U are char- 
acterized by small positive 7-values, revealing that the 
potential drop at the soft contact (e. g., AFM-tip) is 
slightly larger, 8V t > 8V S [cf. Eq. ((2])]. So, even a (very) 
small difference in the interfacial potential drops causes 
a significant polarity dependence of the transition volt- 
age. This suggests that, more than the active molecule 
(which can differ, see Table|l|, the contacts are important 
for the Vt±-asymmetry and calls for a systematic exper- 
imental investigation on the role of the contact groups 
(thiol, amine, etc). 

The 7-values of Table Q] are significantly smaller than 
those estimated via DFT calculations (7 ~ 0.2 — 0.3)^ 
and that of 7 = 0.25 claimed 2 ^ to be appropriate for 
the experiments of Refs. [l] and [IH With the value 
7 = 0.25, Eq. © yields HOMO-offsets e h = 0.87 eV 
(for V t+ = 0.62 V) and e h = 0.97 eV (for V t+ = 0.69 V) 
for the I-F-curves shown in Refs.[l| andUH, respectively. 
The very small difference (much smaller than the experi- 
mental inaccuracies) between the above values and those 
given in Ref. [H (s h « 0.85 eV and 1.0 eV), deduced by 
assuming a Lorcntzian transmission [the key assumption 
underlying Eq. (J3pi] , demonstrates again the high accu- 
racy of the present estimates based on Eq. (|4]). 

While the e^-estimates (eh = 0.62 eV and et = 0.69 eV 
deduced in the present paper versus Eh = 0.85 eV and 
Eh = 1.0 eV of Ref. [29h cannot be directly checked against 
experiment, there is a conclusive test to decide which 
potential division factor (7 ~ 0.07 or 7 = 0.25) is correct. 
Namely, for 7 = 0.25, the transition voltages for negative 
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FIG. 5: The experimental J-F-curves measured in CP-AFM 
setup by Beebe et at^ and Fan et a l 12 ' 30 for anthracene- and 
terphenyl-based junctions, respectively plotted along with the 
theoretical curves computed via Eq. Q with the eh- and 7- 
values of Fable [I] Also shown are the theoretical curves de- 
duced by using the value 7 = 0.25 given in Ref. [Isl . which 
completely disagree with experiments for negative voltages. 
See the main text for details. 



bias are V t - = -2.17V and V t - = -2.42 V, and they 
clearly disagree with the experimental data of Refs. [D 
andKL2L respectively. 

It is worth mentioning that, as visible in Fig. 01 the 
Vt+/Vt- asymmetry very rapidly varies with the division 
potential factor 7. It becomes important even for small 7, 
as revealed by Table |U Although more or less asymmetric 
J-V-characteristics [/(V) 7^ — 7(— V)] are ubiquitous in 
molecular electronics, such a large asymmetry (a factor ~ 
3.5) between the positive and negative transition voltages 
(V t+ = 0.62V and V t+ = 0.69V versus V t - = -2.17V 
and Vt- = —2.42 V, respectively) has not been reported 
so far. 

A curious aspect should still be noted at this point. 
By using the values Eh — 0.87 eV and Eh = 0.97 eV (the 
curves are indistinguishable from the choice Eh = 0.85 eV 
and Eh = 1.0 eV), and 7 = 0.25, and by adjusting the 
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prefactor in Eq. (|4|) to reproduce the experimental linear 
conductance, the agreement with experiments for positive 
biases is also good; see Figs. and b. However, as vis- 
ible there, there is a complete disagreement for negative 
biases. 

Although the DFT inability to solve level alignment 
problems is well known ; 32 ' 33 it has been claime d 18 ' 29 that 
DFT-cstimatcs of the ratio eo/V* could be trusted. The 
DFT-values of 7(~ 0.2 — 0.3) are based just on the fact 
that the ratio e /V t is determined by the value of 7. The 
above analysis reveals a further limitation of the DFT- 
based approaches to molecular transport, raising serious 
doubts on its reliability for TVS studies. 

One should finally note that the above quantity Eh 
is to be understood as the HOMO-offset characteriz- 
ing a molecule linked to two (non-biased) electrodes. 
Therefore, it is not surprising that the present values 
(eh = 0.62 eV and Eh = 0.69 eV for anthracenthiol- and 
terphenylthiol-based junctions) deduced from transport 
data can differ from the ultraviolet photoelectron spec- 
troscopy (UPS) datai Indeed, they do differ: UPS mea- 
surements on anthracene- and terphenyl-thiol on gold 
yielded e^ ps ~ 1.7 cV and 1.8 eV X Partly, the differ- 
ence between Eh and can be attributed to image 
effects i 29 ' 34 but I do not address this issue here. 



V. CONCLUSION 

Although the role of the asymmetric interfacial po- 
tential drops has already been noted in the original 
TVS paper— and considered more recently ; 18 ' 19 ' 29 a di- 
rect quantitative analysis of full i-U-curves measured in 
experiments has not been attempted in previous studies. 
The present paper shows that, by resorting to ambipolar 
TVS, it is possible to determine not only the energy offset 



£0 (at which the original TVS aimed), but also the poten- 
tial division factor 7 without a notably increased effort: 
the only quantities to be determined experimentally are 
the transition voltages Vt± for both bias polarities. For 
quantitative estimates, simple, very accurate analytical 
formulas have been deduced, which have been excellently 
validated against available experimental data. The excel- 
lent agreement found in the present paper gives a strong 
support to the correctness of the Lorentzian transmis- 
sion [the key assumption underlying Eqs. fl3J) and (@J] and 
rules out possible higher-order contributions (~ V 2 ,V 3 , 
etc) to the RHS of Eq. Q| fo r the molecular junctions 
investigated in Refs. Q] andli~2l. 

The present analysis has emphasized the need to in- 
clude both bias polarities in order to obtain correct Eq- 
and 7-estimates. In this context, a specification for 
the discussio n 13 ' 35 of experimental data is helpful. Ac- 
cording to a result of previous work^ the ratio be- 
tween the energy offset and the transition voltage |eo|/Vt 
varies from 0.86 to 2. The present paper reconfirms 
this result with an important specification: the afore- 
mentioned range refers to ratio \£o\/Vt^ m im which cor- 
responds to the transition voltage of the smallest mag- 
nitude V t ,min = min (V t+ , —V t -). The ratio \so\/V t , m ax, 
which corresponds to the transition voltage of the largest 
magnitude Vt^max = max (Vt+, — Vt—), ranges from to 
0.866 [cf. Eq. © and Fig. Si]. 
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